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Abstract: We apply recently developed path integral resummation methods to perturba- 
tive quantum gravity. In particular, we provide supporting evidence that eikonal graviton 
amplitudes factorize into hard and soft parts, and confirm a recent hypothesis that soft 
gravitons are modelled by vacuum expectation values of products of certain Wilson line 
operators, which differ for massless and massive particles. We also investigate terms which 
break this factorization, and find that they are subleading with respect to the eikonal am- 
plitude. The results may help in understanding the connections between gravity and gauge 
theories in more detail, as well as in studying gravitational radiation beyond the eikonal 
approximation. 
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1. Introduction 

Perturbative quantum gravity has been well-studied over many decades, in the form of 
general relativity minimally coupled to matter particles (see [1] for a recent review), as 
well as alternative theories. Although the pitfalls of such an approach are well-known - 
chiefly that quantum general relativity contains non-renormalizable ultraviolet divergences 
- it is nevertheless useful for a number of reasons. Firstly, one may treat perturbative 
quantum gravity as an effective theory, and calculate loop corrections to gravitational 
observables. Secondly, it may well be the case that if a consistent quantum theory of gravity 
is found, it shares features with quantized general relativity. This motivates the study of 
features in this theory which, it may be argued, should be generically shared by alternatives. 
A third reason for studying gravity is that there are intriguing connections between the 
structure of perturbative scattering amplitudes in gravity and non-abelian gauge theories 
(e.g. [2, 3, 4, 5, 6, 7]). A possible application of these ideas is, for example, that of settling 
the question of whether N = 8 supergravity is ultraviolet finite [8, 9, 10, 11, 12]. 

In light of the above remarks (and as stressed recently in [13]), an interesting prop- 
erty of perturbative gravity is the structure of infrared singularities. One might hope, for 
example, that the infrared behaviour of alternative (possibly UV finite) quantum grav- 
ity theories shares at least some of the features that underly quantum GR coupled to 
matter. Furthermore, the structure of IR singularities in renormalizable gauge theories is 
well-understood. IR divergences cancel between real and virtual contributions for suffi- 
ciently inclusive observables [14]. The singularity structure nevertheless becomes relevant 
in that it affects residual large contributions to the perturbative expansion of physical 
cross-sections which must be resummed to all orders in the coupling constant in order to 
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obtain sensible physical predictions. Resummation is by now a highly developed subject 
in both abelian [15] and non-abelian gauge theories. A number of approaches exist, such 
as Feynman diagram methods [16, 17], Wilson line techniques [18, 19], effective field the- 
ories [20, 21, 22, 23, 24, 25], and the recently developed path integral approach of [26]. 
Crucial to all these approaches is the notion of factorization, namely that scattering ampli- 
tudes separate into a hard interaction, which is infrared finite, and a soft function, which 
contains all infrared singularities. These are generated by the emission of soft (zero momen- 
tum) gauge bosons between the external lines of the amplitude, referred to in the literature 
as the eikonal approximation^ . Such factorization properties are important both for resum- 
mation applications, and also for understanding the structure of infrared singularities to 
all orders [33, 34, 35, 36, 37, 38, 39, 40, 41]. 

The same physics also governs the infrared properties of perturbativc gravity ampli- 
tudes, as first analysed in [32]. Recently, Naculich and Schnitzer have reconsidered the 
structure of gravitational IR divergences [13]. The main point of their paper is to argue 
that infrared singularities at all orders are generated by the exponentiation of the one-loop 
divergence i.e. that there are no subleading divergences. Their argument is based upon 
the assumption that gravitational amplitudes factorize in the same way that gauge theory 
amplitudes do, in terms of hard and soft functions. That is, the n-graviton scattering 
amplitude (adopting the notation of [13]) may be written 



where is IR-finite, and Sn collects all IR singularities, generated by soft graviton ex- 
change. They make the further hypothesis that, by analogy with gauge theory, the soft 
function can be expressed as a vacuum expectation value: 



where there is one such operator associated with each external line of the amplitude (with 
momentum pi). Here /i^'^ is the graviton field (we will define this more carefully in what 
follows), and k = 327rG, where G is Newton's constant. The line integral in eq. (1.3) is 
over a straight line contour, parametrized by s. 

The aim of this paper is to examine these hypotheses within the path integral resum- 
mation framework of [26] , which was developed in the context of abelian and non-abelian 

^When amplitudes involve massless external particles, hard coUinear singularities must also be taken 
into account in so-called jet functions [27, 28, 29, 30, 31]. Those are not usually considered in gravity, as 
coUinear singularities cancel after summing over all emitting particles [32]. 

^These operators describe the emission of soft gravitons from hard emitting particles, and should not 
be confused with the parallel transport operator of general relativity, defined in terms of the Christoflel 
symbol, which is also sometimes referred to in Wilson line terms. 
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containing the Wilson line operators^ 
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gauge theory. The main purpose of that paper was to address corrections to the cikonal 
approximation, and the result was a systematic classification of the next-to-eikonal con- 
tributions to scattering amplitudes - that is, those which occur at first subleading order 
in an expansion of the amplitude in the momenta of emitted gauge bosons. The results 
were subsequently confirmed using Feynman diagrammatic methods in [42] , and have also 
been used to classify the structure of (non-abelian) soft gluon corrections in multiparton 
processes, generalising the concept of webs from two parton scattering [43, 44, 45] to the 
multiparton case [46, 47]^. The results of [26] show that the structure of a scattering 
amplitude A subject to soft photon emissions has the schematic form 



Here Aq is the hard interaction amplitude undressed by soft photons, and the exponent 



and next-to-eikonal order respectively. A set of effective Feynman rules for forming these 
diagrams has been given in [26, 42], which generalize the well-known eikonal Feynman 
rules to subleading order in the momentum expansion. The term Arem. in eq. (1.4) is also 
next-to-eikonal order, but does not formally exponentiate. It consists of diagrams in which 
eikonal photons are emitted from an external line and land inside the hard interaction. We 
may call these internal emission diagrams to distinguish them from the external emission 
diagrams which enter the exponent, and their origin makes clear that they can be thought 
of directly as a breaking of the factorization of the amplitude into hard and soft parts. 
Although internal emission contributions do not formally exponentiate, they do have an 
iterative structure to all orders in perturbation theory, which is fixed by gauge invariance 
by a result known as the Low- Burnett- Kroll theorem [57, 58] (see also [59] for a further 
generalization). Examples of internal and external emission contributions are shown in 
figure 1. 

Although the case of abelian gauge theory is described in the preceding paragraph, a 
very similar structure occurs in non-Abelian theories (essentially, the sum over connected 
diagrams in eq. (1.4) is replaced by a sum over so-called webs, as described in [43, 44, 45, 
48, 46, 47]). We will see in this paper that the same general structure is also observed in 
perturbative gravity, assuming that the infrared properties can be analysed independently 
of the uncertain UV completion of the theory. In particular, we will derive directly the 
appropriate form of the Wilson line operator which describes soft graviton emissions, as well 
as the appropriate gravitational generalization of Low's theorem, which shows that internal 
emission contributions in gravity are subleading with respect to the eikonal approximation, 
as is also the case in gauge theory. These results strengthen the hypotheses made in [13], 
and may also shed more light on the correspondence between gauge theory and gravity. 

There are also phenomenological motivations for studying corrections to eikonal grav- 
ity. The latter has been used in a variety of recent applications, in particular the study 

^See also [48] for an alternative viewpoint on multiparton webs. Other work on next-to-eikonal correc- 
tions can be found in [49, 50, 51, 52, 53, 54, 55, 56]. 
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contains connected subdiagrams (which span the external lines) and at cikonal 
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Figure 1: Examples of (a) a generic hard interaction with outgoing particles, which may emit soft 
gauge bosons; (b) an external emission contribution; (c) an internal emission contribution. 

of transplanckian scattering in different gravitational theories. This has a number of po- 
tential uses, such as investigating whether different gravity theories have the same long 
distance behaviour [60], or even what potential collider signatures are in extra dimension 
scenarios [61]. The classification of corrections to the eikonal approximation, as examined 
in this paper, may allow further study of these and related subjects. 

The structure of the paper is as follows. In section 2 we describe the path integral 
framework for scalar particles emitting soft gravitons, adapting the approach used for 
abelian gauge theory in [26], and derive the form of the Wilson line operator for soft gravi- 
tons. In section 3 we discuss how gauge invariance can be used to constrain factorization- 
breaking terms, and discuss how this relates to Low's theorem [57] in QED. In section 4 
we discuss our results before concluding. 

2. Path integral approach to soft graviton amphtudes 

In this section, we will consider the scattering amplitude for L scalar particles, dressed by 
any number of soft gravitons. Although pure multigraviton amplitudes were considered 
in [13], scalar particles will be sufficient for our purposes, due to the fact that eikonal 
gravitons are insensitive to the spin of the particle which emits them [32]. 

The action for general relativity minimally coupled to a (complex) scalar field is given 



here -Se.h. is the Einstein-Hilbert action (which must be suitably gauge-fixed in order to 
define the graviton propagator), and^ 



in d dimensions, where g is the determinant of the metric tensor g^^'^. Perturbation theory 
can then be defined after expanding the metric tensor about the fiat space Minkowski 

''Note tliat we use the metric (-,+,+,+), as in [26]. 



by 



Sgrav = SE.H.lgn + -S„,at [0*, 0, ^''1 , 
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metric rj'^'^ . As explained in e.g. [1], there is an ambiguity in how one performs the weak 
field expansion. Here we adopt the approach of [62], and define 



(2.3) 



This introduces a Jacobian in principle in eq. (2.2), although this may be taken to be one 
in dimensional regularization [1]^. We then define the graviton field h^^'^ via 



(2.4) 



where = 32ttG as in [13]. Note that symmetry of the metric tensor implies h'^'^ = h^^^. 
One could also have chosen, of course, to expand the metric g^^^ directly. However, the 
choice of eq. (2.4) results in simpler expressions for the scalar-graviton vertices. We will 
return to this point later on when discussing next-to-eikonal corrections. 
With the above definition for h^^ one has [62] 



——h" -4- «-2 



a\2 
a) 



d-2 



L(d-2)2 

Furthermore, the inverse of g'^'^ is given by [63] 

The action of eq. (2.2), evaluated to quadratic order in k, is thus 



(2.5) 



(2.6) 



— m^0*^ — 



a\2 



d-2 ' ^(^^-2) 
For what follows it is useful to write this as 

S = - j d'^x(j)*S(f>, 
where the quadratic operator S is given by 



(2.7) 



(2.8) 



+ 



d-2 



K + K 



{d-2) 



(2.9) 



(n.b. we have integrated by parts where necessary, neglecting surface terms). 

The starting point of the path integral approach of [26] is to separate the gauge field 
(here the graviton) into hard and soft modes. That is, in the path integral which defines 
the quantum gravity theory, one may write 



(2.10) 



^Many other field redefinitions arc possible instead of that of eq. (2.3). See [2] for a recent proposal in 
the context of exploring the factorization of gravity amplitudes into gauge theory amplitudes. 
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where and h'i'^ contain only hard and soft modes in momentum space respectively. The 
exact details of this definition (which may proceed e.g. by constructing an explicit hyper- 
surface in the multigraviton momentum space which separates hard and soft parts) need 
not concern us, as discussed in [26] for the QED case. However, we note that this separa- 
tion breaks the full gauge invariance of the gravity theory. Under a gauge transformation, 
the graviton behaves as 

h^^'^ix) h^'^ix) + d^eix) + d'^eix) (2.11) 

where x denotes 4-position. A general gauge function (defined in position space) will 
in general have both soft and hard modes in momentum space, so that gauge transforma- 
tions exist which violate the separation into soft and hard modes introduced in eq. (2.10). 
However, a residual gauge invariance remains, namely that one may transform h^^^ in mo- 
mentum space according to any transformation^ 

Ks^k) h^,^^{k) + k'^a,sik) + k^e,,sik) (2.12) 

in which the gauge function contains only hard or soft modes as required. This sym- 
metry will be sufficient to derive the structure of factorization-breaking corrections from 
internal emission graphs in section 3. 

We now consider the Green's function for the scattering of L scalar particles which, 
after performing the above separation of the graviton field, may be written 

G{pi . . . ,Pl) = / Vhi^'^Hixi, . . . ,XL)e^^^ « ['^n l[{pj\{S - ie)-'\xj). (2.13) 

Here H{x\, . . . ,xl) is the hard interaction, which produces scalar particles at 4-positions 
{xi}, whose external momenta are given by {pi}- This is directly analagous to the hard 
interaction given for the QED case in [26], except for the fact that this will now consist 
of a sum of Feynman diagrams involving hard graviton modes, rather than hard photon 
modes. There are also integrations over the positions {a::^}, which we do not show explicitly 
in eq. (2.13). Associated with each external line is a propagator for a scalar particle 
in a background soft gravitational field. According to the well-known definition of the 
propagator, this is given by the inverse of the quadratic operator S for the scalar field (f), 
which includes the relevant gravitational interactions as shown in eq. (2.9). The propagator 
is sandwiched between states of given initial position and final momentum. The latter 
is due to the fact that Green's functions (and, consequently, scattering amplitudes) are 
usually considered in momentum space. The former arises because, in what follows, we 
will interpret soft graviton scattering in terms of the space-time worldlines of the emitting 
scalars which participate in the hard interaction. 

The next step is to note that the propagator factors in eq. (2.13) can be expressed in 
terms of first quantized path integrals. Such a technique arose some years ago [64, 65, 66], 

®We use the same symbol to denote the position- and momentum-space graviton fields, where the 
argument of the field removes any ambiguity. 
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motivated by string theoretic approaches to field theory amphtudes [67, 68] (see also [69] 
for an example of worldline techniques applied in a resummation context). Here we shall 
quote the result as written in [26], that the propagator factors appearing in eq. (2.13) may 
be written 



{Pj\iS-iey 



where 



1 roo rp(T)=Pj 
— 0^1 ~ \xj) = - dT / PpPxexp — 

2 JQ Jx{0)=Xj 

+i j dt{p ■ X — H{p, x)) 
Jo 



H = -S. 



ip{T) ■ x{T) 



(2.14) 



(2.15) 



Some explanatory comments are in order. Firstly, there is a double path integral over the 
position and momentum trajectories of the particle, which are parametrized in terms of 
a time-like variable t (the upper limit of T will eventually be taken to oo, corresponding 
to the final state). The boundary conditions are that the particle must be produced at 
position Xj (a;(0) = Xj), and end up with final momentum pj {p{T) = pj). Equation (2.14) 
arises as the solution of a Schrodinger equation for a certain evolution operator formed 
out of S, and indeed has the recognizable Feynman path integral representation, where 
the second term in the exponent is the classical action formed out of the "Hamiltonian" 
H. The first term in the exponent arises due to the fact that we are sandwiching the 
propagator between position and momentum states, rather than two position states. Note 
that the path integral over x in eq. (2.14) has a clear physical interpretation in terms of 
a sum over all possible spacetime trajectories of the scalar particle whose final momentum 
is Pj. This will be crucial to isolating the properties of the eikonal approximation (and 
beyond) in what follows. 

From eqs. (2.9) and (2.15), we find that the appropriate Hamiltonian operator for a 
scalar particle coupled to gravity is (dropping the subscript s on the soft graviton field h't^ , 
which we do from now on) 



1 o m^K 
+ -m^ + 



2'" ' 2(d-2) 
which in momentum space (p^ = id^y becomes 

1 



(d-2) 



(2.16) 



H=^-{p'+ m") + ^p.pji^'' + '-^Pu{d^hn 



+ 



2{d - 2) 



(2.17) 



'^Care must be taken with the sign of the momentum operator, which is here chosen to ensure consistency 
with the diagrammatic calculation of appendix B. 
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so that the propagator function of eq. (2.14) becomes 

rp{T)=pj 



{pj\{S-ie) ^\xj) 



1 poo 

'^l 



VpVx exp < —ip{T) ■ X 



,{T) + i f 
Jo 



dt 



p ■ X 



rn^K 



2{d - 2) 



{hi 



a\2 



Pa 



(2.18) 



Xd-2) 

At this point, we can recover the eikonal approximation. Recahing that this corre- 
sponds to the momentum k of any emitted gravitons being completely soft (A; — > 0), this 
means that the emitting scalar particles do not recoil. Thus, they follow the straight line 
classical trajectories 

x(t) = Xj+pjt. (2.19) 

Subeikonal corrections correspond to a systematic expansion about this trajectory. That 
is, for each external line one writes 



X{t) Xj + pjt + x{t), p{t) Pj + p{t), 



(2.20) 



where the boundary conditions for the transformed variables are x{0) = p{T) = 0. One 
then finds [26] 



{Pj\{S-ie) 



-^\xj) = - r dTe-'P^-''i~^^i+"''^^fiT), (2.21) 
2 Jo 



where 

fiT) 



p(r)=o 



VpVx exp < i 



a;(0)=0 



dt 







P-X-'^P^ - '^{Pjf, + P^){pj,y + Pu)h'"' 



K + K 



a\2 



(2.22) 



Xd-2) 

The path integral over p is Gaussian, and may be performed explicitly. First, one writes 
eq. (2.22) as 



J X 



p(T)=0 



VpVx exp 



'a;(0)=0 

where we have defined 



. Jo 



dt 



-p^A^'^p, + Bi'p^ + C 



(2.23) 



IK 



2 



h'^ + K 



(2.24) 



Xd-2) 

and also used the fact that h^'^ = W^. Using standard results for Gaussian integrals gives 

1 



f{T) = y" Par exp i £ dt (^^B''{A-^)^,B'' + C 



(2.25) 
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where we have ignored an overall normalization constant (which is cancelled upon correctly 
normalising the measure for the path integral). Noting that A^'^ as defined in eq. (2.24) 
is just the metric tensor g^'^ of eq. (2.3), its inverse is simply {A^^)^^, = g^y, as given in 
eq. (2.6). Then eq. (2.25) gives (after some tedious algebra) 



/(T) = j Pxexpji j dt 



K in 

— —X tlixuX — Kpjail - 



2 -2 2 



2 {d-2) " 2 \{d-2y d-2 



(2.26) 



One may relate all this back to the eikonal scattering amplitude as follows. After 
combining eq. (2.26) into eq. (2.21) and substituting the resulting propagator factors into 
the Green's function of eq. (2.13), one finds (after truncating the external lines with free 
scalar propagators as required by the LSZ formula) that the scattering amplitude for L 
scalar particles has the form (see the derivation in [26]) 

/oo 
Vh^^'H{xi, . . . , Xi)e^^^-H-['^''"'l n e-^^'^-^V.(oo). (2.27) 

That is, each external line is associated with a factor /j(oo), given by eq. (2.26) with 
r ^ oo, which still implicitly includes a path integral over all possible trajectories of 
particle j. We will shortly analyse this result further, but first note that eq. (2.27) has the 
form of a generating functional for a quantum field theory, i.e. for the soft graviton field. 
The external line factors act as source terms in this theory which, as can be seen from 
eq. (2.26), generate soft graviton emission vertices localized along the external lines. The 
Feynman diagrams generated by eq. (2.27) are thus subdiagrams in the full theory, which 
span the external lines. Given that connected diagrams exponentiate in quantum field 
theory, we can immediately conclude that a large class of soft graviton corrections (namely 
those consisting of connected subdiagrams formed from the Feynman rules derivable from 
eq. (2.26)) exponentiate. This is a generalization of the original result in [32], which is 
based solely on the eikonal contributions. The argument given here for the exponentiation 
of soft graviton corrections based on the textbook exponentiation properties of connected 
diagrams in quantum field theory is essentially exactly the same as the argument used for 
QED in [26]8. 

Let us now physically interpret the various terms in the external line factor of eq. (2.26) 
in more detail. First, note that the Feynman rules generated by each term cannot be 
immediately read off, due to the fact that for each external line the path integral over 



^Strictly speaking, exponentiation of connected diagrams holds only subject to the condition that the 
source terms for the gauge field are mutually commuting. This is true for gravity and abelian gauge theories, 
but is not the case in non-Abelian gauge theories, where a more sophisticated treatment of exponentiation 
is then required [26, 46, 48]. 
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trajectories x{t) has yet to be carried out. However, as shown in [26], this can be done 
by systematically expanding about the classical trajectory. The leading term (as already 
remarked above) gives the eikonal approximation, in which the particles do not recoil. The 
first subleading corrections give the next-to-eikonal corrections, and further corrections are 
next-to-next-to-eikonal etc. Put another way, an expansion about the classical trajectories 
of the hard emitting particles (in position space) is entirely equivalent to an expansion in the 
momenta of the emitted gravitons (a momentum space expansion). We may evaluate the 
eikonal contribution in eq. (2.26) by setting x{t) = p(t) = 0, and also neglecting terms which 
are quadratic in h^^ (these lead to two-graviton emission vertices, which are suppressed 
with respect to the eikonal limit, due to having one less eikonal denominator compared with 
two separate graviton emissions). One may also neglect any terms containing derivatives 
of the graviton field, as these (in momentum space) will be suppressed by powers of the 
graviton momentum. The result is 



/f'^(oo) = j Pxexp ij^ 



dt 



(2.28) 



Carrying out the path integral over x gives unity assuming the appropriate normalization, 
and thus the eikonal approximation for the scattering amplitude is 

7=1 L ^0 



m 



2{d-2) 



(2.29) 



In fact, up to next-to-eikonal corrections, one may set the initial positions of the external 
lines to xj = [26], so that eq. (2.29) simplifies further to 



/oo r /■ 
p/,M-giSB.H.[h^1-rTexp i 
.1=1 L Jo 



dt ( -^Pji.Pj.h'''' (2.30) 



m 



2{d - 2) 



(2.31) 



where we have taken the hard interaction outside the path integral over the soft gauge 
field. One thus sees that the amplitude has factorized explicitly into a hard and a soft 
part, where the soft gravitons are described by the Wilson line factors 



= exp 



oo 

I I dt 





I"-"-""" - 2(1^"-'"'"' 



(2.32) 



The generating functional for the soft graviton field theory indeed generates vacuum expec- 
tation values of these operators, as expressed in eq. (1.2). Let us now discuss the physics 
of these operators in more detail. 

Firstly, the first term in the exponent in eq. (2.32) agrees, up to a normalization 
factor, with the hypothesized Wilson line operator of [13] (which dealt with pure graviton 
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scattering, and thus massless external particles). Note, however, that the Feynman rule 
for the soft graviton generated by this term is derived from 

-'if = / (0 (5^"-(^) . (2.33) 

where h^'^{k) is the momentum space graviton field. The momentum space eikonal Feyn- 
man rule is thus 

2p-k' 

whose normalization agrees with the result quoted in [13]. Thus, we maintain that eq. (2.32) 
is correctly normalized. 

Secondly, note that there is an additional contribution in the case in which the external 
particles are massive, as has been considered here, which appears as the second term in 
the exponent of eq. (2.32). Wilson line operators in conventional gauge theories are the 
same for massless and massive particles, due to the fact that the mass does not source the 
gauge field. Here, however, the mass plays the role of a gravitational charge, and thus we 
expect such a dependence. Furthermore, such a term can be motivated on general Lorentz 
invariance grounds, given that the exponent of the Wilson line operator can only depend 
on the properties of the emitting particle. That is, the prefactor which one contracts with 
the graviton field can only depend upon and m, and the only two combinations having 
two Lorentz indices and the correct mass dimension are 

which indeed both occur. 

Thirdly, it is interesting to note the behaviour of eq. (2.32) under rescalings of external 
line momentum. In QED, the Wilson line operator has the form 



^QED = exp 



ie J dtpf.A'^ 



(2.34) 



where e is the electromagnetic charge. This leads to the momentum-space eikonal Feynman 
rule 

^ (2.35) 
p - k 

for emission of a soft photon of momentum k, and this latter expression is manifestly invari- 
ant under the transformation p — t- Xp, corresponding to the invariance of eq. (2.34) under 
reparamctrizations of the Wilson line. By contrast, the eikonal Feynman rules generated 
from eq.(2.32) do not have this property. Rather, they can be written in the form 



K 

2 



, Pu \ , 'Htiu m 



(2.36) 



The quantities in the round brackets are both manifestly scale invariant under rescalings 
of the external momentum. These are then multiplied respectively by factors of p^ and 
m, namely by gravitational charges. Thus, the eikonal Feynman rules in gravity have the 
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form of a scale-invariant (in p) quantity multiplied by a charge, which is exactly what one 
expects by analogy with abelian and non-abelian gauge theories. 

The curious reader may wonder about the presence of poles in d — 2 in the Wilson 
line operators. As discussed in e.g. [1], perturbation theory becomes singular in d = 2 
dimensions, owing to the fact that the Einstein-Hilbert action is then a topological invariant 
(i.e. there is no non-trivial dynamics). Whether or not such poles occur in the Wilson line 
operators is dependent on the choice of weak field approximation. If one expands the 
metric itself rather than absorbing y/—g, the {d — 2) poles instead occur in the graviton 
propagator, such that one can never get rid of them, as expected from their physical origin. 

Above, we have confirmed the factorization property of soft graviton amplitudes in 
the eikonal approximation, in terms of vacuum expectation values of appropriate Wilson 
line operators. Some comments are in order regarding the plethora of additional terms in 
eq. (2.26). As explained in [26], the path integral over x in eq. (2.26) can be systematically 
expanded about the leading eikonal term. This was done for the case of abelian and non- 
abelian gauge theory by defining p = An, where A — >■ oo isolates the eikonal term, and 
= or = 1 for massless and massive particles respectively. The path integral could 
then be done perturbatively in A"-*^ (also taking care to expand the argument of each gauge 
boson field), leading to a set of next-to-eikonal Feynman rules. There were two types 
of next-to-eikonal vertex: those involving a single photon or gluon emission, and those 
involving two gauge bosons emitted from the same position on the external line. These 
vertices were then confirmed by an explicit diagrammatic treatment in [42]. 

A similar procedure is possible in the present case. Note, however, that the simple 
scaling pj = Xnj does not solely isolate the eikonal term as A ^ oo. Instead, eq. (2.26) 
becomes (as A oo) 



lim f{T) = JvxexpS^i J dt 



(d - 2)2 d-2 



(2.37) 



The first two terms are the eikonal contribution of eq. (2.32), whereas the last two term 
are two-graviton vertices, which have no analogue in the case of abelian or non-abelian 
gauge theory. They are clearly suppressed in the eikonal limit, as such vertices have one 
less eikonal denominator than two separate eikonal emission vertices. The above remarks 
imply that the A scaling of [26] is not quite correct as a systematic method for expanding 
about the eikonal approximation^, although does indeed lead to correct results for the 
gauge theories studied in that paper. The extra terms arise here essentially because of the 
weak field expansion, which allows the possibility to construct terms which are quadratic in 
external momenta or masses, and quadratic (or higher) in the graviton field. Such terms are 
absent in QED, where no weak field expansion is needed. In any case, it is straightforward 
to isolate next-to-eikonal contributions from power-counting, even if one does not use the 



^We will see in what follows, however, that the third term in eq. (2.37) is cancelled upon performing the 
path integral over the trajectory x. 
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A-expansion (which amounts to httle more than a convenient book-keeping device where 
apphcable). To calculate the next-to-eikonal graviton Feynman rules, one must perform 
the path integral perturbatively as explained in appendix B of [26]. Here we perform this 
calculation in appendix A, with the result that, up to next-to-eikonal order, the external 
line factor of eq. (2.26) may be written (taking T oo) 



/(oo) = exp 



2 p-k 



+ 



4 d-2{p-kf'^^\ 



A p-k 



V/j.i' -|- ^ 



2 {d-2)p-k A{p-kf 



+ 



, PaPl3P{^ilu) , , , , 



PtiPuPaPl3 {k ■ I) 



PaP0P{^iK) 



+ 



p ■ Ip ■ (k + l) 

V^ii^ P(a 1 13) 



16 \p ■ kp ■ (k + l) 



p ■ Ip ■ {k + 1) 



8 p-kp-lp - (k + l) m{d-2) \p- kp - (k + l) 

Vfii^ PaPi3 [k ■ I) 



4 2 

m K 



+ 



8{d-2) \p-kp-lp- {k + l) 



8(d -2Yp-kp-lp-{k + l) 4p-{k + l) \ id-2)^ d-2 



+ {fiu -H- a/3) 



(2.38) 



Here we have written the exponent explicitly in momentum space, where we use the ar- 
gument of the graviton field to identify its Fourier-transform. Prom eq. (2.38) one can 
easily read off the effective Feynman rules for graviton emission up to NE order. There are 
eikonal and NE one-graviton vertices given by 



K 1 



2p-k 



PixPu + 



m 

d-2 



(2.39) 



and 



Pi/iku) 



+ 



k^ 



p-k ■ (p-k^ V"^"^ d-2- 

respectively, where the eikonal result for massless particles is already well-known [32]. 
There is also an effective two-graviton vertex, given by 



m 



(2.40) 




K 



IQp - kp-lp- {k + l) 

2m2 



Tfl \ 



d-2 



rif,vPaP^{k ■ I) + {l^f o aj3) 



+Arr?{p -k){p- 1) 



id - If 



1 



d-2 



+ 2p^p^PaPl3{k ■ I) + 



{rifMlSrii'a + riu/3'nij.a) 



2m^ 



{d-2)- 



:VtMiyVai3{k ■ I) 
(2.41) 



One expects such a vertex for two reasons. Firstly, whilst successive graviton emissions are 
completely decoupled in the eikonal approximation, this is not expected to hold beyond the 
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eikonal approximation, such that one expects two-graviton correlations. Secondly, there is 
an exact two-graviton vertex in the theory (as discussed in appendix B), which forms part 
of the result of eq. (2.41). It is useful to cross-check the above results using diagrammatic 
methods, as has been carried out for the case of QED and QCD in appendix B of [26]. We 
present such a calculation for the present case in appendix B of this paper. 

The above results imply that a subset of next-to-eikonal corrections exponentiates at 
NE order in perturbative quantum gravity. Namely, connected external emission graphs 
formed from effective eikonal and NE Feynman rules, where each graph contains at most 
one NE vertex. Some further comments are in order regarding the results we have obtained. 

To start with, one may question how useful the NE Feynman rules are in practice. 
Eikonal Feynman rules tell us about infrared singularities, which may be expected to be 
somewhat universal in alternative field theories of gravity. Ncxt-to-eikonal corrections are 
non-singular, and thus potentially tell us less about generic infrared features of gravity 
theories. Furthermore, the fact that perturbative GR coupled to matter is not UV renor- 
malizable suggests that an expansion in the momenta of emitted gravitons may cease to 
be meaningful at some point. Thus, one should be wary of taking any NE corrections 
too seriously. What matters in the above analysis is merely that in the soft limit, we 
have shown that the eikonal approximation leads to Wilson line operators as hypothesized 
in [13], where factorizablc corrections arc genuinely sublcading in that they do not produce 
IR singularities. We have also clarified the situation when external particles have non-zero 
mass. 

Another important point is that unlike in abelian and non-abelian gauge theory, the NE 
Feynman rules are not unique. They depend on the nature of the weak field approximation 
used to define the graviton field (e.g. whether one expands g^'^ or g^'^). This is not the case 
for the eikonal terms in four dimensions. We have verified by explicit calculation (which 
we do not consider worth recording here) that the Wilson line operator which occurs upon 
expanding g'^'^ rather than g'^'^ is 



Comparing this with eq. (2.34), we see that the two Wilson line operators are equal for 
d = 4. They are also automatically equal for massless external particles. The next-to- 
eikonal terms, however, would differ between the two calculations, due to extra terms 
originating from the expansion of ^/—g in the kinetic energy term for the scalar field. 
Such differences presumably cancel at the level of full diagrams (as the choice of weak 
field approximation shuffles contributions between different vertices), but it is still true 
that a level of ambiguity occurs in thinking about NE Feynman rules that is absent in 
conventional gauge theories, as a weak field expansion is not necessary in the latter. 

A further comment is in order regarding the weak field approximation. In the above 
analysis, wc have performed two expansions. Firstly, we have separated the graviton field 
into hard and soft modes, and used this as a basis for an expansion in the momenta of 
emitted gravitons. Secondly, we have expanded the graviton field itself in the gravitational 
coupling constant k. In the latter expansion, we kept terms only up to 0{k^). How 




(2.42) 
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(a) (b) 

Figure 2: Example external emission graphs which exponentiate at next-to-eikonal order, where 
• represents a next-to-eikonal vertex, the doubled gluon line a graviton, and all other vertices are 
assumed to be eikonal. 

can we be sure that by including higher order terms, we do not induce further (next-to) 
eikonal corrections? To see that this is not the case, it suffices to note that any power 
of K is accompanied by a further power of the graviton field h^'^ , whose Lorcntz indices 
must be appropriately contracted with other factors of /i^*^, pj, or in cq. (2.26). 
This always results in a multigraviton vertex, which is suppressed with respect to the 
appropriate number of eikonal 1-graviton vertices due to the fact that it lacks at least one 
eikonal denominator. The fact that all vertices at 0{k"') generate n-graviton vertices also 
shows that for next-to-eikonal corrections it is sufficient to expand only up to O(k^), as we 
have done above. This is itself interesting, as it shows that the momentum expansion and 
weak field expansions are partially correlated in a well-defined sense (although there is, of 
course, a tower of subleading momentum terms at any given order in k). 

It follows from the above discussion that the NE Feynman rules will contain 1-graviton 
and 2-graviton vertices. It then follows (as for the eikonal corrections) that a large class 
of next-to-eikonal graviton corrections exponentiates, namely those graphs containing con- 
nected subdiagrams involving one next-to-eikonal Feynman rule. Examples are shown in 
figure 2. Thus, the exponent of the scattering amplitude for L-particle scattering dressed 
by soft gravitons has the same schematic form as the QED case of eq. (1.4). As in that 
case, this is not the whole story. One may also get additional contributions from internal 
emission diagrams such as that shown in figure 1(c), which explicitly break the factorization 
of the scattering amplitude into hard and soft parts. These are known to be subleading 
in the soft limit in conventional gauge theories [57, 58, 59]. It is likely, but not perhaps 
obvious, that the same is true for gravity. This is the subject of the following section. 

3. Factorization breaking terms 

In the previous section, we applied the path integral resummation method of [26] to per- 
turbative general relativity minimally coupled to a complex scalar field. We found that 
the amplitude factorizes into hard and soft parts, where the latter is decribed by a vacuum 
expectation value of Wilson line operators as in eq. (2.32). Explicit calculation yields dia- 
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grams formed out of cikonal Feynman rules, which then exponentiate. Corrections to the 
exponent are subleading, and consist of diagrams containing one next-to-eikonal Feynman 
rule, with all other vertices eikonal. The aim of this section is to examine what hap- 
pens when this factorization is broken, and to argue that, as in QED, such corrections are 
next-to-eikonal, thus do not affect the structure of IR singularities. 

First let us briefly recap what happens in QED. We denote by H^^{pi, ...pL;k) the sub- 
amplitude for production of a soft photon emerging from the hard interaction H{pi, . . . p^) 
with momentum k and with Lorentz index fi. Then gauge invariance may be applied to 
show that H'^ is related to the hard interaction itself, up to next-to-eikonal order. More 
specifically, this relation has the following form: 



where qj is the electric charge of particle j. This equation tells us that factorization 
breaking terms are subleading with respect to the eikonal approximation. Furthermore, 
their structure is completely fixed up to next-to-eikonal order. Such contributions do not 
formally exponentiate (as do external emission graphs), but they do nevertheless have an 
iterative structure to all orders in perturbation theory, as given by eq. (3.1). This equation 
is an important ingredient in Low's theorem [57, 58, 59], which relates the vertex function 
in gauge theory to the vertex function with an additional emission of a gauge boson. 

The above result suggests that we should attempt a similar strategy in quantum gravity. 
That is, one may impose gauge invariance on the factorized form for the soft graviton 
scattering amplitude, and check that this leads to extra contributions which are suppressed 
with respect to the eikonal limit. As discussed in the previous section, the separation of the 
graviton field into hard and soft modes has broken the full gauge invariance of the theory. 
Nevertheless, the residual gauge symmetry under transformations of the type shown in 
eq. (2.12) will be sufficient to derive a gravitational analogue of eq. (3.1). 

Our starting point is the factorized form of the soft graviton amplitude of eq. (2.29), 
where each external line factor is reduced to a Wilson line operator. Under a gauge trans- 
formation (eq. (2.11)), the Wilson operators transform as 



Requiring that eq. (2.29) be gauge invariant then implies that the hard interaction trans- 
forms as 




(3.1) 



(3.2) 




(3.3) 



Expanding both sides to first order in h^^ and ^'^ gives 

Hixi,...,XL) + J d''xH^,{xi,...XL;x) [h''''ix) + d''eix) + d''eix) 
= H{xi,...xl) + J (fxHf,^{xi,...,XL\x)h^"'{x) 



H{xi, . . . XL)^PjnS{x - Xj 



(3.4) 



where all hard functions are evaluated with the untransformcd graviton field hM'^ , and 
H^i,{xi, . . . ,xl;x) is the subamplitudc for emission of a graviton from the hard interaction 
at position x (this is the analogue of H^{xi, . . . ,xl;x) in the QED case, whose Fourier 
transform appears on the left-hand side of eq. (3.1)). Equation (3.4) simplifies to 



J Ah^^xu . . . xi;x) la^Cix) + 8T(x)l 

=»/A;/(..,...x.)i:p,,*(x-x,)ew, 

j 

and integrating by parts on the left-hand side gives 

- J dV^r^^, = m J d'^xH{xi,...XL)Y,P3Ax-X3)i'"{x)^ 

3 

where H^j^i, = H^j^y + Hy^. Due to the fact that is arbitrary, one may write 



(3.5) 



(3.6) 



-S^iJ^i, = iK^pj^H{xi, . ..xl)5{x - Xj). 



(3.7) 



In momentum space this becomes 



(3. 



and expanding this to 0{k) gives 



, ^ d 



H{pi,...,pl) 



(3.9) 



The zeroth order term on the right-hand side cancels due to momentum conservation 

J]p,-^ = 0, (3.10) 



and one finally finds 



d 

H^,u = -i^^Pjf^QpijHipi, ■ ■ ■ ,Pl)- 

3 ^ 



(3.11) 
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This is the gravitational equivalent of the QED result expressed by eq. (3.1). It forms part 
of a gravitational version of Low's theorem, and describes the result of internal emission 
diagrams analagous to that shown in figure 1(c), in which a soft graviton emerges from 
the hard interaction and lands on an external line. Clearly such diagrams are next-to- 
eikonal from eq. (3.11). To see this, imagine starting with an external emission graph, 
and taking one of the soft gravitons off an external line and into the hard interaction. 
Such a replacement replaces an eikonal Feynman rule with the contribution of eq. (3.11), 
and the resulting diagram thus has one less eikonal denominator of form p ■ k, therefore 
is subleading with respect to the eikonal approximation. It is worth noting that one may 
also derive eq. (3.11) from a traditional diagrammatic treatment - we present an example 
calculation in appendix C. 

There is also an additional contribution at NE order, corresponding to the fact that 
the external lines do not start at the origin, but have non-zero initial 4-positions Xj, as 
taken into account by the factors e~^i'^i in eq. (2.27). As in [26], one may evaluate 
this contribution by writing the eikonal one-graviton source term for an external line of 
momentum p as 



2 J (2^ 



Ax-k 



2 J (2^ 



d'^k hi"" 



Yp-k 



m 



Pt^ Pu + -J—jPt^^j i^ + ix-h), 

(3.12) 

where we have expanded to NE order on the right-hand side. Combining the second term 
on the right-hand side with the rest of eq. (3.11), the scattering amplitude is given by 



d\ 
? 



A{p,,...,pl) = I me^^-.M (^~i:Pj,^H{p,...,p,)^ I ^ 

+ Jdxi... dxiH{xu ...Xn)lY.j (^(^^^- • . . k 



^Pj^^Pju + ^Vf.. )\e ^^•^'^ 



N 



(3.13) 



where we have explicitly reinstated the integrals over the initial positions of the external 
lines, and f{0,Pj; h) is the external line factor starting at the origin for line j. Performing 
the position integrals, one obtains 



N 



'^Hipi,...,pL)}llf{0,pj;h). 



(3.14) 



One sees that the result of the non-zero initial positions of the external lines is to add an 
extra NE term involving momentum derivatives of the hard interaction. A similar result 
occurs in QED and QCD, and it is useful to physically interpret this extra contribution. 



-18- 



In abclian and non-abelian gauge theory, massless particles give rise to collinear as well 
as soft singularities, and one may write the scattering amplitude in the factorized form 
(analagous to eq. (1.1)) 

L 

A = H-S-^Jj. (3.15) 

Here H and S are the hard and soft functions, where the latter collects all singularities due 
to the emitted gauge boson momenta satisfying A; ^ 0. Furthermore, Jj is a jet function, 
and contains all singularities arising from hard emissions which are collinear to the outgoing 
particle of momentum Pj^^- Del Duca has shown [59] that in such cases there are extra 
contributions to the Low-Burnctt-KroU theorem due to emissions of soft gluons from inside 
the jet functions. This is indeed the physical meaning of the additional term in eq. (3.14). 
One would not ordinarily write jet functions in perturbative gravity, due to the fact that 
hard collinear singularities cancel after summing over all external particle lines and using 
momentum conservation [32]. However, if one is interested in the structure of ncxt-to- 
eikonal corrections, there are such terms which do indeed come from emissions within jet 
functions, and which do not cancel after summing over parton legs (due presumably to the 
fact that the additional term in eq. (3.14) is quadratic in each external momentum). It 
may thus be useful to consider jet functions in gravity after all. 

In this section, we have investigated the structure of terms which break the factoriza- 
tion of soft graviton amplitudes into hard and soft parts, and found contributions whose 
physical meaning is that of emission of soft gravitons from inside the hard and jet func- 
tions. All such contributions arc sTibleading with respect to the eikonal approximation, 
thus verifying the factorization of graviton ampitudes at eikonal order. Furthermore, the 
structure of soft graviton amplitudes up to NE order has the same schematic form as in 
QED or QCD, namely that of eq. (1.4). External emission graphs (consisting of connected 
subdiagrams formed out of eikonal and next-to-eikonal Feynman rules) formally exponenti- 
ate, whereas internal emission graphs have an iterative structure and do not formally enter 
the exponent. 

4. Discussion 

In this paper, we have studied the recent hypothesis made in [13], that soft graviton ampli- 
tudes factorize into hard and soft functions, as do gauge theory amplitudes. Our analysis 
relied on path integral resummation techniques developed in the context of abelian and 
non-abelian gauge theory [26], which allow an efficient classification of the structure of 
scattering amplitudes including subleading corrections with respect to the eikonal approx- 
imation. We find that gravity amplitudes have the same schematic structure as in gauge 
theory, as written here in eq. (1.4). That is, eikonal graphs exponentiate (as already 
known [32]), and next-to-eikonal graphs can be separated into external and internal emis- 
sion contributions. External emission graphs are connected subdiagrams formed out of 

^°Note there is an overlap between S and Jj , namely when singularities become soft and collinear. This 
is usually dealt with by dividing either the soft or jet functions by the eikonal jet functions evaluated in 
the soft limit, which removes any double counting. 
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effective eikonal and ncxt-to-eikonal Feynman rules, wliose form we derived in this paper. 
Internal emission graphs consist of soft gravitons which are emitted from inside the hard 
interaction, or jet functions associated with each external line. External emission graphs 
formally exponentiate, whereas internal emission graphs have an iterative structure to all 
orders in perturbation theory. Importantly, internal emission graphs are subleading with 
respect to the eikonal approximation. 

Note that, as in the cases of QED and QCD considered in [26], one may quibble 
whether it is really valid to separate internal and external emission contributions as far as 
exponentiation is concerned, or even what is meant by next-to-eikonal exponentiation at 
all. After all, upon performing the exponent of the NE external emission contributions, 
any term involving products of two or more NE terms is NNE order or higher, such that 
eq. (1.4) is formally equivalent to 



in which formula external and internal emission contributions are on the same footing. 
Nevertheless, it is true that the former genuinely exponentiate, whereas the latter have an 
iterative structure and thus must be calculated differently. Furthermore, either eq. (1.4) or 
eq. (4.1) allow NE contributions to be calculated to all orders in perturbation theory, even 
if this is really a consequence of the exponentiation of eikonal terms, as shown explicitly 
in eq. (4.1). 

Aside from the structure of NE effects, an important question regarding these results 
is: have we formally proved the factorization of graviton amplitudes (albeit those involving 
only scalar external lines) into hard and soft parts? After all, the path integral technique 
we have used starts off by assuming that we may separate the gauge field into hard and 
soft modes, and that the resulting exponential structure captures all singularities. It may 
perhaps be prudent to regard the results we have obtained as being supporting evidence for 
the factorization property of amplitudes, rather than a rigorous proof. Nevertheless, it is 
very strong evidence, for a number of reasons. Firstly, we have analysed corrections which 
explicitly break the factorization property (i.e. internal emission graphs), and found these 
to be subleading in the momenta of the emitted gluons. We showed this using both the 
path integral framework, and also using diagrammatic methods (although this is essentially 
the same calculation, as it is gauge invariance that ultimately fixes these contributions). 
Secondly, we find that the structure of eikonal and next-to-eikonal contributions is exactly 
analagous to that of QED and QCD. This suggests that the qualitative structures one finds 
in the infrared limits of both gauge and gravity theories are somewhat generic. In gauge 
theories, the factorization property (together with the structure of the NE Feynman rules) 
is indeed explicitly corroborated using diagram-based methods, which gives us confidence 
that the same would be true for gravity. A more formal proof of the factorization of gravity 
amplitudes may perhaps proceed using a power-counting analysis based upon the Landau 
equations associated with the singularities of Feynman diagrams, as has been carried out 
in gauge theories (see e.g. [70] for a pedagogical review). 




(4.1) 
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We are also able to confirm the hypothesis of [13] that the soft graviton part of the 
scattering amplitude is described by a vacuum expectation value of certain Wilson line 
operators, with one such operator associated with each external line. This in itself is not 
so surprising: in the eikonal approximation external particles do not recoil, and thus can 
only be changed by a phase. If this phase is to have the correct gauge transformation 
properties to slot into an amplitude, then some sort of Wilson line operator is inevitable. 
In gravity, the Wilson line operator can only depend upon p^Pu and m^rj^u due to Lorentz 
invariance. It is therefore reassuring to have verified that this is indeed the case. The fact 
that massive external particles give an additional contribution to the Wilson line operator 
is also expected, given that and m play the role of charges in gravity. The exponent 
should consist of charges multiplying factors which are invariant under rescalings of the 
external momenta, which is indeed the case. 

In this paper we have considered scalar particles emitting soft gravitons. One could 
also consider the more phenomenologically relevant examples of fermions, gauge bosons 
or hard gravitons emitting gravitational radiation. This should make no difference to 
the eikonal approximation, which is known to be insensitive to the spin of the emitting 
particles [32]. Beyond the eikonal approximation, and by analogy with the QED case 
examined in [26, 42] one expects NE corrections which are insensitive to the spin of the 
emitting particle, together with additional contributions which couple graviton emissions 
to the spin of the emitter. For fermions, this would need an appropriate curved space 
formulation of the second order fermion action [71], which was used in [26] to analyse 
fermions in the path integral approach to soft gauge boson emission. 

The results presented here provide further insight into the relationships between gauge 
theories and perturbative quantum gravity. They may also be useful in phenomenological 
applications of gravitational physics e.g. extending results obtained in the eikonal approxi- 
mation for transplanckian scattering [61, 60]. The investigation of such higher order effects 
in both gravity and gauge theories is ongoing. 
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A. NE Feynman rules from the path integral approach 

The aim of this appendix is to explicitly carry out the path integral of eq. (2.26), perturba- 

tivcly about the classical trajectory. We will keep terms up to and including next-to-eikonal 
corrections, for which a weak field expansion to 0{k^) is sufficient. The calculation closely 
follows that of the QED case in appendix B of [26] . 
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Wc begin by noting that the first term in the exponent of eq. (2.26) defines the 
quadratic operator (after taking T — oo) 

where we have integrated by parts on the right-hand side. The propagator for the 
coordinate is given by the inverse of this operator, which is [26] 

G{t,t') = imm{t,t')rjf,^,. (A.2) 

This in turn leads to the following correlators, which we will need in what follows: 



{x^{t)x^{t')) = G{t,t') = imin{t,t'X 
{x^{t)x,{t')) = ^^Ml = io^t' - t)v^, (A.3) 

{x,{t)xS')) = = iS{t - t')r,,., 

where Q{t — t') is the Heaviside function. Note that the second and third correlators vanish 
at equal times {t = t'). To see this, one may use finite difference formulae consistent with 
the space-time discretization assumed in the path integral derivation of [26] to write 

{x,it)x.{t)) = hm M)Mt + ^)--^m ^ 0, (A.4) 
e->0+ e 

where we have used the first correlator of eq. (A.3) for the final equality. Likewise, one 
may write 

{x^{t)x,{t)) = -{x^,{t)x,{t)) (A.5) 
(using integration by parts), which gives 

{x,{t)xm = - lim Mt)lMt + 2e) - x.jt + e) + x^t)]) ^ ^ 
e->-0+ e 

We now consider the external line factor of eq. (2.26), where for brevity in what follows 
we write pj — ^ p. Furthermore, we may use the symmetry of h'^'^ to rewrite eq. (2.26) in 
the T ^ oo limit as 



/(oo) = J Vx exp |i J dt 



x^ K . , „„ in ,^ , ,,„, . 



+PuPam^,+P^Paripu)h^''h^^ - ^p^p.h^^" - YP{-d^)h'"' - y^^^M-^'^" 



(A.7) 

In this expression, we have anticipated that some of the terms in eq. (2.26) will not lead 
to NE contributions after carrying out the path integral, and so have already discarded 
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them (e.g. terms with both powers of x and more than one graviton tensor). Ignoring 
terms which depend exphcitly on m?' for the moment, one may expand the arguments of 
the graviton tensors (each of which is evaluated at 4-position pH + x'*) where relevant to 
give 



/(oo) = j Pxexp^i j 



dt 



—Pf^Puidadrh'^n^'^X- - '-^Pi.d^)h'''' - '-^Pi^.d^^id.hnx" 



IK 



+-rP(u%))3h^"idah' 



(A.8) 
(A.9) 



where all graviton tensors are now evaluated on the classical trajectory. Also, we have 
neglected terms of form x'^x'^ and x'^x'^, which would vanish in the following calculation 

due to the equal time commutators of eqs. (A. 5) and (A. 6). It is now convenient to 
transform to momentum space (in order to ultimately be able to read off the momentum 
space Feynman rules). Writing 



h^^ix) 
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ik-p t 



(A.IO) 



(where as in the rest of the paper we use the same symbol for the Fourier transformed 
graviton field), the graviton contributions to the exponent in eq. (A.9) can be written 

{ / (^^'"'^^^^'^ (^-'^PiuX^,) + Ih^X^) - '^Pf.pr, - "^P^Pu Xa 
+ -:PnPuk k XaXr + -7P(uki^) + -rP{uk^)k X„ 
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(A.11) 



In this expression, the terms with no x powers are 



I e 
'0 



ip-k t 



+Pi^PaVPiJ.+PiJ,PaVpu) 



PiJiPv + -^P{vk^) - -re''^'^^P{yrii^)liK + — e*^'''* {PuP/dVai^ + PiJ.Pi3r]au 



(A.12) 



where to shorten the notation (and from now on) we omit the integrals over momenta 
and the Fourier transformed graviton fields, which can easily be reinstated based on the 
number of Lorentz indices in each term. The various terms in cq. (A.12) can immediately 
be interpreted in terms of Feynman rules, as the remaining path integral over x gives unity, 
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assuming the appropriate normalization. Terms with one power of generate vertices for 
the coordinate x^, which we label as follows: 



'0 



[A) 



^ • ^ 1 • J cr 7 7 (J 

"^Piu^H) ~^~r'^(iJ,Xu) ~~^PtiPu k Xcr +-rP(i/"^/i)^ 



(A.13) 



In carrying out the path integral over x to NE order, one must construct all tree level 
graphs containing two such vertices. Again omitting the graviton fields and momentum 
integrals, we start by considering the (A)-(A) combination, which tree level graph gives 



1 / in 



dt 



/ dt'e'^-'''+P-'''\\p.x^{t)+p^x,{t)]\ppXa{t') +paxp{m, (A.14) 
Jo 



where the initial factor of 1/2 is a symmetry factor. Inserting the correlator of eq. (A.3), 
eq. (A.14) becomes 



J dt e*P'('=+')* \pi,pi3 rjan + Pi^PlB rjau + PuPa + P/xPa Vpv] ■ (A.15) 

Note that this precisely cancels the contribution from the final term in eq. (A. 12), which 
thus does not result in a contribution to the effective Feynman rules at NE order. 
Next, one has the (A)-(B) combination, which gives 

(~y) ^'^■'^^'"'^'^(b. i^it) + P^ (*)] [fcai/3(t') + kpx^it')]) 

= y dt e'P'^^^^>* \py ka r],j,i3 + ka ?7i//3 + Piykj3 r]^a + Pii k/3 flva] ■ (A.16) 

Exploiting the fact that there is an implicit factor of h^'^ {k)h"'^ {I) in eq. (A. 12), we may 
further rewrite eq. (A.16) as 



IK 



2 i-oo 

I dt e*P-('=+^)* [p^ ka rjij.^ + Pu. ka + Puk^i ViJ.a + Pn kn r]ya + {k^ I)] . (A.17) 

^0 

By the same token, the third term in eq. (A. 12) can be written 



^^2 roo 

~ le" / e*^''(*^+')* \pi, ka rif,i3 + Pfj, ka r]^^ + Pv kp r]^a + Pix kjj rj^a + {n,i',k ^ a, P, I)] 

(A.18) 

The terms shown in eq. (A.18) precisely cancel eq. (A.17), and thus there is again no 
contribution to the effective Feynman rules. 
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Next up is the combination (A)-(C), which gives 



poo poo 

Jo Jo 



4 p • kp • {k + I) 

r P<xPpP{fiK) , PaPpP{tiK) , , „ ■ 

H , , + [iiv ^ a/3) 



ie> [p- kp- {k + l) p-lp-{k + l) 



(A.19) 



where have symmetrized over graviton indices in the final Une. This is indeed a NE con- 
tribution, and thus contributes to the effective two-graviton vertex. 

By a similar procedure, one may evaluate the contributions from all other combinations 

of the vertices in cq. (A. 13). The reader may verify that (A)-(D), (B)-(B), (B)- (C), (B)- 
(D), (C)-(D) and (D)-(D) do not give any contribution up to NE order. For (C)-(C), one 
has 

Pi^PuPo^P^k'^f dt dt'x^{t)xr{t')e'^-'''+P-'''^ 
Jo Jo 

^■^P^P.Po.Pp k-l dt dt' ram{t, ty^-^t+rp-it') 

o Jo Jo 

PlxPuPaPpk-l 



2 V2 



8 p - kp-lp - {k + l)' 



(A.20) 



This is already symmetric under (iiv a/3). 

One must also consider the term in eq. (A.ll) with two powers of x. This contributes 
to a loop graph at NE order, which gives a contribution (including a symmetry factor of 
1/2) 



1 foo 

—PfiPuk'^k^ / dtXa{t)xr{t)e'P-''' 

= -^Pt.Puk' r dtte'P-^* 
_ Kp^puk'^ 



[p ■ ky 



(A.21) 



This is already symmetric in /i and u. Furthermore, the fact that this term is quadratic in 
x means that the symmetry factor of 1/2 does not enter the resulting Feynman rule. 
The uncancelled terms with no x powers from eq. (A. 12) give 



IK 
-^Pi^Py 



IK 



JO 



dte'P-''' = - 



2 p-k' 
k P(h k^) 
A p-k ' 



(A.22) 
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Above we have ignored dependent terms, and it is straightforward to repeat the 
above analysis for these. Expanding arguments of graviton tensors where necessary and 
keeping only relevant terms, the mass-related terms in eq. (A. 7) give 



dt 



4 d-2 



2 \{d-2)-' 2{d-2) j'^ \^ 



(A.23) 



which in momentum space becomes 



[{I 



d'^k 
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{27r)d 
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im K 



2 d-2 



Y d-2 



2^2 „ {r]ij,fsVua+Viyl3Vtia) 



{d-2) 



d-2 



hi"'{k)h'^''{l) } . 
(A.24) 



The terms with no powers of x give (again omitting the graviton fields and momentum 
integrals for brevity) 



dt e 



ip-kt 



2 d 



f 

Jo 



2 2 

im K 



2 d-2p-k' 



2 2 

m K 



4p- {k + I) 



2 

2ri^u riai3 

(d-2)2 {d-2) 



_(d-2)2 2(d-2) 



1 



The term in eq. (A.24) which has one power of x gives a vertex for x'* 



m 



2 d 



— I e ^ Tjjuj/ k X(j . 
— ^ Jo 



(A.25) 



(A.26) 



(A.27) 



Labelling this by (E), one must consider all tree level graphs in which (E) is contracted 
with one of the vertices (A) to (D) of eq. (A. 13). First of ah, (E)-(A) gives 



m K 
~2 d-2 



) poo 

dt dt'{\paXp{t') +P0Xa]Xa{t)e 

Jo 



i(p-kt+p-lt' 



m 



A{d-2)p-kp- {k + l) 



+ 



+ {i^v <H- aP) 



(A.28) 



16(d -2) [p- kp- {k + l) p-lp-{k + l) 

The reader may check that (E)-(B) and (E)-(D) do not give NE contributions. The 
combination (E)-(C) gives 



2 d-2 
_ W? ■q^jLvPaPp k ■ I 



oo />oo 

dt / dt'x^Xre'^P-'^'+P-^''^ 
Jo 



4 d-2p-kp-lp- {k + l) 

rifMuPaP/B k ■ I 



8 d-2 [p-kp-lp- {k + l) 



+ {iiv -H- aP) 



(A.29) 
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p-k 



Figure 3: Graph used to check the one graviton emission vertex. 



Finally one has (E)-(E) which gives 

Jo Jq 

(A.30) 



2 V 2 d-2^ 



8{d-2)p-kp-lp- {k + l)' 

There is also a term in eq. (A.24) which has two powers of x. This gives a loop contribution 
to the path integral 



2 4 d-2 



Vi^.k'^k-J^ dte'^-''\x,{t)xr{t)) = —j-^rj^,j^, (A.31) 



where again the symmetry factor of 1/2 does not enter the Feynman rule due to the fact 
that this term is quadratic in x. 

This completes the calculation of all contributions to the x path integral up to NE 
order. Collecting the results of eqs. (A.19, A.20, A.21, A.22, A.25, A.26, A.28, A.29, A.30, 

A. 31), one indeed finds that the external line factor is given by eq. (2.38). 

B. Diagrammatic check of the NE Feynman rules 

In this appendix, we cross-check the effective Feynman rules for graviton emission up to 
next-to-eikonal order, using a diagrammatic approach. Our analysis is similar to that 
presented in appendix B of [26], and rests upon the fact that the effective Feynman rules 
must reproduce the known NE limit of an exact calculation which includes up to two 
graviton emissions. 

We begin by considering figure 3, which contains a single graviton emission from an 
external line of momentum p. The exact Feynman rule for the graviton emission is given in 
eq. (C.l), where pi and p2 are the momenta to the left- and right-hand side of the vertex. 
Given that this assumes all momenta are outgoing^^, one must set p\ = —{p — k) and 
P2 = p- Combined with the propagator for the internal graviton line, this gives a total 
contribution 

-2 {p-kY + m- + 0^2"^^^) ■ ^^-'^ 



^See also footnote 6 on p. 7. 
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p-k-l 



p-k-l p-l p-k-l p-k 

(a) (b) (c) 

Figure 4: Graphs used to check the two graviton emission vertex. 



Expanding this up to NE order gives 



+ 



+ 



+ 



m 



2 p-k A p-k 2 d-2p-k 4 {p-ky 4 d-2'^''{p-k)'^' 



(B.2) 



This agrees with the results of eqs. (2.39) and (2.40). 

To check the two-graviton vertex, we consider the graphs shown in figure 4. Firstly, 
there is the graph of figure 4(a), consisting of two separate single graviton emissions. This 
gives a contribution 



T 



(p-fc-0(^(p-/)^) 2m^ ri^y 
{p-k-lf + 'm? d-2{p-k-lf + 'm? 



' {p-l){aPp) 2w? r]a(s 

{p-lf + m? d-2{p- If + m? 

(B.3) ■ 



In checking the two graviton vertex, one may combine this graph with its counterpart 
in which the gluon momenta are permuted, with a factor 1/2 to correct for the double 
counting. The sum should then separate into a part which factorizes into products of 
lower order emissions, and a non-factor izable remainder, which should agree with the two- 
graviton vertex obtained using the path integral approach [26, 42]. The permuted graph 
is shown in figure 4(b), and its contribution can be easily obtained from eq. (B.3) by 
interchanging k and I. 

There are two types of contribution in expanding the weighted sum of these graphs 
to NE order. Firstly, there are terms which arise from expanding to 0{k) or 0{l) in the 
numerator. Secondly, there are terms which result from expanding the denominators about 
their eikonal limit i.e. 



1 



1 



{p — kY + rn^ 2p - k 



1 + 



A:2 



2p-k 



(B.4) 



where the ellipsis denotes terms beyond the NE approximation. Considering first the 
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numerator contributions, expanding eq. (B.3) and combining this with figure 4(b) gives 

\P^i Pv Pa Pli (k + ^■)ifiPu) Pa Pl3 ni^ {-q^i, Pa Pl3 + Va/S Pfj. Pv) 



p- (k + l) 2p- 



+ 



+ 



+ 



d-2 



2{d-2) p-ik + l) {d-2Yp-{k + l) 



p-ik + l) 

1 1 

+ 



p ■ k p ■ I 



k{^iPu) Pa Pp 



{aPp)PtiPv 



{aPp) PtiPu 



m 



klJiPv) Pa Pli 

2p-{k + l)p-l 

VaP kp.Pi') 



2p- {k + l)p-k 2p-{k + l)p-l 2p- {k + l)p-k 2{d - 2) p ■ Ip ■ {k + I) 



2{d-2)p-kp-{k + l) 2{d-2)p-kp-{k + l) 2{d-2)p-lp- {k + l) 

The square-bracketed terms factorize upon apphcation of the eikonal identity 

1 1 
+ 



p ■ {k + I) \p ■ k p-lj p-kp-V 



(B.5) 



(B.6) 



to give products of lower order effective Feynman rules. Left are the non-factorizable 
contributions 



c2 r / 



{lxPv)PaPp 



16 [p- (A; + /)p- Z p-{k + l)p-k d - 2p ■ {k + l)p ■ I d - 2p ■ {k + l)p ■ k 



+{fj,i' O a/3) 



(B.7) 



These thus contribute to the effective two-graviton vertex. 

Now considering the denominator terms in eq. (B.3) (together with the corresponding 
result with k I), the NE terms from the m^-independent terms are 



K^P^iPuPaPp / {k + l) 



1 1 

+ 



8 p-{k + l) \2p-{k + l)\p-k p-l 

which can be rearranged to give 

k!^ P^iPuPaPp l'^ 



+ 



k' 



+ 



+ 



I^^PfiPu PaPp 

16' 



+ 



+ 



2{p-k)^ 2{p-l)^ 

PfiPuPaPpk ■ I 



p-l {p ■ k)^ 16 p-k {p ■ l)^ 8 p ■ kp ■ Ip ■ {k + I) 



(B.8) 



(B.9) 



The first two terms have factorizcd into products of lower order effective Feynman rules. 
The final term, however, is non-factorizable and hence contributes to the effective two- 
graviton vertex. 

Now considering the m^-dependent denominator contributions, these are given by 
(keeping only NE terms) 



m? {ri^iuPaPp + rjapPuPu) , rri 



_d-2 p-{k + l)p-l 



+ 



nixv Va/3 



{d-2)^p- {k + l)p-l 



+ 



{k + 1)' 
2p-{k + l) ' 2p-l 



+ 



(A; o Z) } . 



Rearranging and keeping non-factorizable terms only, one finds 



K^m'^ {rifj.^PaPl3 + 'nal3PfiPu)k ■ I 

8{d-2) p-kp-lp- {k + l) 
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1 4 

Km 



Vfiiy Vapk -l 
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(B.ll) 
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Figure 5: (a) Example vertex function involving two outgoing particles; (b) Series of diagrams 
involving all possible placements of an extra graviton. 



Finally, there is the graph of figure 4(c), which is given by the exact two-graviton 
vertex [1] 

r 2 ] 

riu0 + VfiP Vua - J—^^f^'^ VaiS ■ (B-12) 



2{d-2) ^ 

combined with an eikonal propagator, to give 



2 2 

Km 



1 



A{d-2)p-{k + l) 



d-2 



(B.13) 



Adding this to the non-factorizable contributions from eqs. (B.7, B.9, B.ll), one finds 
a total effective two-graviton vertex at NE order that is equal to that presented in eq (2.41) 
(obtained from the path integral approach). One could of course have used the diagram- 
matic method of this appendix to derive the NE vertices in the first place, without relying 
on the path integral approach. However, there is then no proof that the effective Feynman 
rules hold to all orders in the perturbation expansion i.e. that they are genuine Feynman 
rules^^. This latter point is made manifest by the path integral method. 



C. Diagrammatic derivation of eq. (3.11) for internal emissions 

In this appendix, we give an example of how eq. (3.11) (which expresses how the subam- 
plitude for soft graviton emission from within the hard interaction can be related to the 
hard interaction with no additional emission) can be derived from a traditional Feynman 
diagram calculation. 

We start with the vertex function shown in figure 5(a) with two external scalar parti- 
cles, which in figure 5(b) is dressed by an additional graviton emission in all possible places. 
For simplicity, we deal with the case of massless external particles. The vertex with no 
emission can then only depend on the invariant pi ■ p2, and we denote this by T{pi,p2). 
The full Feynman rule for emission of a graviton from a scalar is given by [1] 



d-2 



(C.l) 



^An all order diagrammatic proof of the NE Feynman rules in QED and QCD has been given in [42]. 
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in the weak field approximation of eq. (2.4). For m = the vertex function including an 
extra graviton emission (from figure 5(b)) is then 



r 



K 

4 



-2pi ■ k 
P2ix{P2 - k)y + (p2 - k)^,P2v 



■T[{pi - k) ■P2] 



-2p2 ■ k 



rbi • (P2 - k)] 



I T-iint. 



(C.2) 



where FJf*' represents the contribution from the internal emission graph. It will be conve- 
nient to symmetrize on both sides, defining F^,^ = F^j, -t-Fj,^ (and likewise for F™* ). Then 
eq. (C.2) can be rewritten 



2 2 



PiniPi -k)„ + {pi- k)^pu 
-2pi • k 



■T[{pi - k) ■P2] 
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Expanding this to first order in k gives 



If 

2 2 



^fr-p..^-^)+^fF-p,.^-^ 

-Pi - k \ opi-P2j -P2 -k \ dpi 



dv 
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2pi ■ k 2p2 ■ k 



^ pint. 
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where F = T{pi •P2), and we have used the notation pa^^ky) = Pifj,ki, +Pipkfj,. One may now 
use the gravitational Ward identity [74] 



fc'^f = 



(C.5) 



(which follows from general coordinate invariance) as well as the on-shell condition for the 
emitted graviton {k!^ = 0) to find 



\ jLMf int. I ^ 

2 '^'^ + 2 



Ply [V -p2-k 



dV 
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P2u - Pi ■ k 
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(C.6) 



The zeroth order term vanishes from momentum conservation (pi +P2 = 0). Also, one may 
identify 

dV dV dV dV 



P2fM 



Opi ' P2 Op\ Opi • P2 Op'2 



(C.7) 



so that eq. (C.6) becomes 



T^int. _ _^ 



dT 



(C.^ 



This is indeed a special case of eq. (3.11) for internal emission graphs. 
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